Dimensional crossover of the fundamental-measure functional for parallel hard cubes 
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We present a regularization of the recently proposed fundamental-measure functional for a mix- 
ture of parallel hard cubes. The regularized functional is shown to have right dimensional crossovers 
to any smaller dimension, thus allowing to use it to study highly inhomogeneous phases (such as 
the solid phase) . Furthermore, it is shown how the functional of the slightly more general model of 
parallel hard parallelepipeds can be obtained using the zero-dimensional functional as a generating 
functional. The multicomponent version of the latter system is also given, and it is suggested how 
to reformulate it as a restricted-orientation model for liquid crystals. Finally, the method is further 
extended to build a functional for a mixture of parallel hard cylinders. 
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Although nowadays the success of density functional 
theories to describe nonuniform fluids and their phase 
transitions is out of question, it is generally acknowledged 
that these theories are based on much to empiric assump- 
tions to be considered a systematic tool to study liquids 
[jlj. Nearly all of the available functionals propose a de- 
pendence on the density through nonlocal weights, whose 
choice rely upon approximations around the uniform fluid 
combined with heuristic arguments. Accordingly, they 
always need a good knowledge of the bulk structure and 
thermodynamics of the uniform fluid, which thus be- 
comes an input of the functional. There is, however, 
a kind of functionals which do not follow this general 
framework and are constructed only from geometric fea- 
tures of the particles; these are the so called fundamental 
measure functionals (FMFs) Two important fea- 

tures distinguish FMFs from "classical" functionals: (1) 
the thermodynamics and structure of the uniform fluid 
(and even that of the nonuniform fluid) can be derived 
from them, instead of requiring it as an input, and (2) 
they achieve a dimensional reduction of the system by us- 
ing strongly inhomogeneous density profiles HH . Actu- 
ally these two properties are interlaced, as it will become 
clear in what follows. 

The first FMF was proposed for a fluid of hard spheres 
(HS) J|||, and it was obtained by assuming it to de- 
pend on a small set of geometrically-weighted densities 
(how many and which weights is dictated by the low 
density limit of the free-energy functional). The precise 
functional dependence still required an ad hoc scaled- 
particle assumption, but in turn, it yielded the Percus- 
Yevick (PY) free energy and direct correlation function 
(DCF), and allowed to compute further structural prop- 
erties (such as the third order DCF [^J). The func- 
tional admits two-dimensional (2D) as well as a one- 
dimensional (ID) ^ extensions, the latter being the well- 
known exact functional. 

There is a feature the exact HS functional exhibits: a 



right dimensional crossover, or, in other words, the fact 
that the functional for D — 1 dimensions must come out 
as a result of the one for D dimensions when the den- 
sity profile is a delta function along one of the coordi- 
nates. But this turns out to be a much to stringent re- 
quirement to be fulfilled by an approximate functional. 
Needless to say that functionals with this property will 
provide good descriptions of strongly inhomogeneous sys- 
tems (e.g. fluids confined by walls). However, it would 
not be of great importance to describe bulk phases were 
it not because there is a special dimensional crossover 
any functional willing to exhibit a sensible solid phase 
should verify: the reduction to a 0-dimensional (0D) sys- 
tem HQ, i.e. a cavity able to contain at most one particle 
(which clearly mimics the situation of a particle in a crys- 
tal lattice). None of the "classical" functionals verify any 
dimensional crossover |9| , but the FMF for HS does some 
of them |^] ; however, when trying to reduce from 3D to 
0D, some divergences appear that cause this functional 
not to be able to stabilize a solid phase. Although an 
heuristic modification of the functional eliminates those 
unphysical divergencies || , it is very remarkable that the 
D-dimensional HS FMF, for any dimension, can be ob- 
tained from the exact 0D functional by requiring it to 
have a right dimensional crossover from D to quasi-0 di- 
mensions H (according to Ref. 0, a quasi-OD system 
is provided by a density function consisting of a sum 
of deltas placed at points such that HS simultaneously 
placed at them overlap) . The results turn out to be non- 
local functionals of the one-particle weighted densities, 
which can be further approximated by local ones without 
loosing their main dimensional crossovers. The moral of 
this is that the 0D functional seems to be the only thing 
one needs in order to make a FMF for a given system. 

Recently, a FMF has been developed for a system of 
parallel hard cubes (PHC) [Q. Although such a system is 
rather unrealistic, it is far more suitable than the HS sys- 
tem to study the demixing transition of mixtures [pULM, 
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and so, it deserves consideration. On the other hand, on 
a fundamental viewpoint, having FMFs for systems other 
than HS may help to clarify the basic structure of such 
functionals in the aim of extending them to more gen- 
eral particle shapes (arbitrarily orientable, anisotropic 
particles, for instance, so as to study liquid-crystalline 
phases). In this letter we will show how the functional 
proposed in |jj can be regularized to have right dimen- 
sional crossovers, and this will make clear a simple rela- 
tion between the FMF for PHC in different dimensions 
which will permit to derive them very easily from the OD 
functional. 

Any FMF expresses the excess (over the ideal) free en- 
ergy of the system in terms of the local density p(r) as 



F cx [p] =k B T J dv^ D \{n a (v)}), 
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where D is the dimension and the n a (r) are a set of suit- 
ably chosen geometrically-weighted densities. According 
to Ref. the free energy density for the system of 



(3D) PHC, is given by <f>^ = + $ 2 3) + $ 
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Here n a = p ® uj a and n Q = p ® w Q1 where <g> denotes a 
convolution. The weights ui a and w Q are given by 
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with t x = Q(a/2 - |a;|), ( x = (\/2)8(a/2 - \x\) (a being 
the edge- length of the cubes), and u = (1,1,1). No- 
tice the different normalization in the definition of the 
weighted densities with respect to that of Ref. [Q . 

As mentioned in [Q , the 2D functional can be obtained 
in a similar way; in this case the weights are given by 



wo(r) = (x( y , w 2 (r) 
wi(r) = {T x ( y ,(xTy) : wi(r) 
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with u = (1, 1), and <J>< 2 ) = $f J + $ 2 2) , with 
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Of course, the ID functional is the same as that of 
HS (i.e. the exact one) because in this case both, HS 
and PHC, reduce to segments on a line. Therefore, 



-noln(l — ni), where = £ x and 
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^i(r) = t x . 

Let us now consider the ID to OD dimensional 
crossover. To this purpose let us define the class of quasi- 
OD densities as the set of functions with support in the 



interval (—a/2, a/2) (the fact that this interval is cen- 
tered around zero is completely irrelevant and it is only 
assumed for simplicity). The densities p(x) = 6(x) and 
P( x ) = J2iViK x - x i)> witn -ct/ 2 < Xi < a/2, r]i > 0, 
and J^iVi = " < 1 (those considered in |^|), belong to 
this class, but smooth functions with support within that 
interval do it as well. Then it is straightforward to check 
that ITl 
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where s x = signs. Now, if we consider that 
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where || = r\ + (1 — r])\n(l — if), we can check 

that F cx = J dx$ w = $ (0) (ni(0)), which reduces to the 
excess free energy of the OD system for any density of the 
quasi-OD class. 

To check the 2D to ID and 2D to OD dimensional 
crossovers, let us define the quasi- ID density class as the 
set of 2D functions with support in the band —a/2 < 
y < a/2, and let us redefine the quasi-OD density class as 
the subset of the quasi-lD class whose members vanish 
for x (—a/2, a/2) (i.e. 2D functions with support in a 
square centered at the origin). Let us consider now the 
2D functional (||), and let us put rix = (ni x ,ni y ); then 
ni — ri! • rix = 2n% , x n\ y . On the other hand, as in the 
previous case, 
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and then 
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which is the right 2D to ID dimensional crossover. Fur- 
thermore, as we have discussed previously, the functional 
( |l0| ) has the right dimensional crossover to OD for any 
quasi-OD density. 

In both cases, ID and 2D, the right dimensional 
crossover arises as a consequence of the fact that for 
quasi-OD densities the functionals can be obtained as 
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The dimensional crossovers then follow from the direct 
integration of the derivatives. It can then be inferred 
that the D-dimensional functional with the right dimen- 
sional crossovers will, for quasi-OD densities, be expressed 
as 
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but if we compute the 3D case, and take into account 
that 
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we find that <&\ and $2 are the same as those in Eq. 
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where the notation v • v • v = v\ + Vy + v 3 has been in- 
troduced. This new expression differs from that of Eq. 
(^J) not only in the numerical coefficients, but also in the 
appearance of a new term, n 2 • n 2 • n 2 , which has no 
equivalent in HS At first this may seem a spuri- 

ous term, because it is not even rotationally invariant; 
on second thoughts, there is no reason why it should be, 
since the system itself (PHC) is not rotationally invariant 
either. Such a symmetry is replaced by a discrete group 
of symmetries which, restricted to the fact that the vec- 
tor densities ni and n 2 must have positive components, 
translates into a simple axes-exchange symmetry. The 
new term obviously has this symmetry. The important 
thing is that while this new 3D functional has been built 

to have right dimensional crossovers, the old one does not 

( 3^ 

even fulfill the 3D to OD one (<I> 3 develops delta cube di- 
vergences when the density approaches a delta function, 
much in the same way as the old version of the HS FMF 
did j|]). It can be checked that the new functional sta- 
bilizes a solid phase, for packing fractions above r\ 0.3, 
by following the standard procedure of parametrizing the 
density function by a sum of gaussians placed at the lat- 
tice sites [Q , whereas the divergences in the old functional 
give rise to a free energy which decreases monotonously 
as the gaussians shrink. 

But to be correct, the new 3D functional should pro- 
vide the exact second and third virial coefficients in the 
density expansion of the DCF, for this is the requirement 
under which FMF are built @,||- It is very remarkable 
that, when applied to the uniform fluid, the new func- 
tional provides exactly the same expressions for both the 



free energy and the DCF as found in Ref. Q (bear in 
mind the different normalization of the weighted densi- 
ties!). Therefore, the new term introduces nothing new 
for the uniform fluid, but it drastically changes the be- 
havior of the nonuniform fluid. 

Formula (|l^) provides a simple method to obtain the 
FMF for a D-dimensional PHC fluid; although the for- 
mula only holds true for quasi-OD densities, once the 
functional form has been derived and the scalar, vec- 
tor, and tensor weighted densities have been identified, 
it can also be applied to arbitrary density functions [this 
is precisely the procedure we have followed to obtain the 
3D functional (|l5j)]. At this level, Eq. ( p"3| ) may seem 
nothing more than a fortunate trick to obtain the func- 
tional; however, a slight manipulation of the equation 
will provide further insight on the actual meaning of this 
method. 

First of all, let us consider the slightly more general 
problem of a fluid of D-dimensional parallel hard paral- 
lelepipeds (PHP) of edge-lengths crW, . . . , along the 
Xi, . . . ,xd axes. If we define rj(r) = p ® uje>(y), with 
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then, the OD free-energy density functional for this sys- 
tem will be given by 



4°^p = ^) + [l-^(r)]ln[l-7 7 (r)]. 
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Formula ([p]) still yields the D-dimensional functional for 
quasi-OD densities, but if we make use of the identity 
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Eq. ( |l3| ) transforms into 
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which is valid not only for quasi-OD densities, but for 
completely general density functions (because it does not 
involve derivatives with respect to the coordinates). 

Equation (^) has the important consequence that 
all the thermodynamics of the inhomogeneous D- 
dimensional fluid is contained in the OD functional. This 
result reveals the same fundamental idea recently formu- 
lated for HS M. In the latter, however, the derivation is 
not as simple as for the PHP system. An attempt to find 
an equivalent to Eq. (18) for HS, where the OD functional 
acts as a generating functional, has just been carried out 
p"2| . In this work, the D-dimensional FMF for HS, in 
the version of Ref. in which all weighted functions 
are scalar, is given by the action on the OD functional 
of a differential operator with respect to the sphere ra- 
dius. The result is nevertheless not as satisfactory as for 
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the PHP system (for instance, the 3D DCF generates 
an unphysical delta function at contact). On the other 
hand, while Eq. ( |l8| ) emerges here in a very natural way, 
the differential operator in [jl2j is introduced heuristically 
with some free constants which are fitted afterwards by 
imposing certain constraints (e.g. to match a given equa- 
tion of state). On the contrary, the functional for 2D 
agrees with the simplified version of the one obtained in 
H , so that the method seems promising and, in view of 
what happens for PHP, it deserves further consideration. 

Generalizing (|l8|) to mixtures is straightforward; it 
simply amounts to writing 
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where denotes the edge-length along the Xj coordi- 
nate of species i. Of course, now 77(1-) = J2i Pi ^w-^'fr), 
where is given by Eq. ( |l6| ) with the a^s corre- 

sponding to species i. The functional ( |l9| ) has a partic- 
ularly important application to study liquid-crystalline 
phases, for if we consider the six possible orientations 
of a parallelepiped as belonging to six different species, 
the functional ( |r| ) can be understood as representing 
a restricted-orientation, one-species system of hard par- 
allelepipeds. Besides that, the functional may provide 
some clues to understand the structure of a FMF for 
freely-orientable, general anisotropic particles. 

The derivation we have made here suggests further ex- 
tensions of the theory. For instance, consider a system 
of hard parallel cylinders oriented along the z axis. It 
is clear that the 3D to 2D dimensional crossover, elimi- 
nating the z coordinate, will transform this system into 
a system of hard disks. Therefore, if <j>^g is a FMF for 



hard disks (e.g. those of Refs. [p|JT2|]), the correspond- 
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ing FMF for the cylinders can be derived from $^3, the 

(2) 

same functional as <&^jg , but with the weights multiplied 
by O (L/2 — \z\) (L being the length of the cylinders) 
and the 2D density of disks replaced by a 3D density of 
cylinders. The derivation will simply be 



$(3) = 
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This functional can immediately be generalized to a mix- 
ture of parallel hard cylinders as in (|l^) , and this provides 
a tool to study analytically the influence of polydispersity 
in the phase diagram of liquid crystals jl3| . 

To summarize, we have shown how the FMF for PHC 
presented in Q can be regularized to have a right dimen- 
sional crossover to any smaller dimension. In the deriva- 
tion we have found the remarkable result that the FMF 



for PHP, in any dimension D, can be obtained by simply 
deriving with respect to every edge-length the 0D func- 
tional, thus prompting the idea of this functional being 
a kind of generating functional (an idea already explored 
for HS |5 12 1, but which for the present system reaches 
its clearest and simplest expression). In passing we have 
obtained the FMF for a mixture of PHP, and suggested 
how it can be exploited as a restricted-orientation model 
for liquid crystals. Finally a further extension of the 
method permits to find a FMF for a mixture of parallel 
hard cylinders, thus opening the possibility of studying 
polydispersity in liquid crystals. 

We are indebted to Y. Rosenfeld for calling our at- 
tention on this problem; to P. Tarazona for illuminating 
discussions and for keeping us informed of his progress in 
studying HS; to J. A. White for sending us their preprint, 
and to B. Mulder for suggesting the extension of the PHP 
functional to a restricted-orientation model. One of us 
(J.A.C.) owes much to a discussion with L. Araujo which 
turned out to be crucial for the development of this work. 
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